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Executive  Summary. 


The  conceptual  design  of  a  novel  electromagnetic  launch  system  has  been  developed. 
While  the  concept  can  be  scaled  in  both  mass  and  launch  speed,  the  performed  study  assumes  a 
100-kg  projectile  leaving  the  launcher  at  7,000  m/sec.  The  projectile  is  propelled  and  levitated 
by  electromagnetic  forces  inside  of  a  low  pressure  tube  filled  with  helium  at  about  10  torr.  A 
fijlly  controlled  acceleration  profile,  which  can  be  adjusted  to  meet  the  requirements  of  various 
missions,  is  used  to  reach  5000  gee  in  a  500-meter  long  launch  tube. 

With  an  individual  launch  cycle  of  approx.  0.5  sec,  several  launches  per  hour  are  possible 
depending  on  the  type  of  energy  storage  system.  Depending  on  range  and  altitude  of  the 
desired  mission,  the  launch  tube  should  be  inclined  to  angles  between  25  and  60  degrees.  To 
achieve  the  desired  launch  angle  with  a  support  structure  easy  to  implement,  the  launch  tube 
can  be  partially  built  underground.  The  system  is  powerful  enough  to  reach  low  earth  orbit,  but 
could  also  be  scaled  to  launch  projectiles  over  distances  of  a  few  hundred  miles  in  kinetic  kill 
applications.  The  proposed  technology  is  highly  reliable  and  responsive:  preparation  for  launch 
(the  time  needed  to  load  the  projectile  into  the  launch  tube,  cool  the  superconducting  coils  on 
the  projectile  and  charge  up  the  energy  storage  system  for  launch)  should  be  in  the  order  of  a 
few  hours  (pre-cooling  the  superconducting  projectiles  would  significantly  increase 
responsiveness  by  enabling  launch  within  less  than  one  hour).  The  system’s  reliability  and 
survivability  is  expected  to  be  very  high,  since  the  projectile  has  no  mechanical  contact  with 
the  launch  track  during  the  acceleration  cycle,  and  the  heat  dissipation  during  launch  is 
minimal. 

The  breakthrough  concept  of  the  proposed  system  is  the  invention  of  a  pulsed 
synchronous  motor  that  generates  a  constant  magnetic  flux  zone  moving  with  the  projectile  in  a 
traveling  wave.  Thrust  and  levitation  is  provided  by  current  pulses  supplied  to  the  track  coils 
adjusted  in  time  (leading  and  trailing  edge,  pulse  duration  and  pulse  frequency)  so  that  no 
magnetic  flux  change  is  seen  by  the  projectile  coils  during  launch.  The  cylindrical  projectile  is 
equipped  with  superconducting  coils  operated  in  persistent  mode,  charged  prior  to  launch. 
Constant  flux  translates  to  no  voltage  or  additional  currents  being  induced,  which  prevents  the 
projectile  coils  from  quenching  or  overheating.  The  projectile  coils,  built  in  a  cable-in-conduit 
assembly  as  in  pulsed  fusion  magnets,  can  carry  persistent  currents  of  up  to  50  kA  under  severe 
operating  conditions  in  pulsed  magnetic  fields  with  no  need  for  induction  or  external  power 
supply.  The  proposed  constant  flux  synchronous  motor  (CFSM)  is  self-centering  and  avoids 
mechanical  contact  between  the  track  and  the  projectile  -  the  electromagnetic  forces  and 
torques  are  sufficient  to  compensate  for  the  aerodynamic  disturbance  forces  and  moments 
during  launch.  The  projectile's  trajectory  is  further  stabilized  by  inducing  spinning  motion  prior 
to  launch. 

A  low  pressure  helium  atmosphere  inside  of  the  launch  tube  is  proposed,  which  nearly 
eliminates  both  aerodynamic  and  heat  transfer  effects  (forces,  torques  and  heating)  due  to  high 
speed  gas  flow  around  the  projectile,  with  large  safety  margins.  The  proposed  technology 
therefore  removes  several  of  the  technical  problems  that  have  made  hypersonic 
electromagnetic  launch  unfeasible,  in  particular  the  inductive  heating  of  the  projectile  when  a 
contact-free  armature  is  desired. 

Injection  to  the  atmosphere  at  ground  level  at  hypersonic  speeds  creates  the  need  of  a 
careful  study  of  ablation  and  heat  management  in  the  projectile.  Our  preliminary  analysis  and 
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existing  studies  indicate  the  feasibility  of  the  proposed  mission;  however,  further  studies  are 
needed  to  folly  qualify  ground-level  launching. 

During  the  proposed  Phase-2  of  this  project,  the  key  components  of  the  system  will  be 
qualified.  This  includes  the  operation  of  the  Constant-Flux-Synchronous  Motor  (CSFM),  the 
superconductor  for  the  projectile  coils  and  the  ablation  and  thermal  management  of  the 
projectile  during  atmospheric  coasting.  The  feasibility  and  efficiency  of  the  CFSM  concept  will 
be  experimentally  demonstrated  on  a  linear  test  rig.  The  required  superconducting  coils  will  be 
developed  and  tested  for  mechanical  and  quench  reliability  under  simulated  operational 
conditions.  Ablation  and  heat  dissipation  in  the  projectile  will  be  tested  using  a  high-intensity 
electron  beam  that  mimics  the  heat  density  induced  by  atmospheric  coasting  -  although  not  all 
aspects  of  the  aerodynamic  erosion  and  heating  during  atmospheric  flight  can  be  simulated, 
important  insight  will  be  gained. 

The  proposed  electromagnetic  launch  system,  based  on  the  constant-flux  synchronous 
motor  concept,  is  a  breakthrough  technology  that  eliminates  several  of  the  inherent  limiting 
factors  of  conventional  electromagnetic  launchers.  The  main  advantages  of  this  concept  are: 

•  The  CFSM  principle  provides  contact-free  propulsion  and  self-centering  suspension 
forces  with  sufficient  stiffness  and  thrust  force. 

•  No  brushes  or  projectile  power  supplies  are  needed,  eliminating  complex  contact 
phenomena  at  high  speeds,  its  inherent  limitation  of  power  rate,  and  undesirable  arcing 
effects  at  the  end  of  the  launcher. 

•  The  use  of  high-current,  persistent-mode  superconducting  coils  in  the  projectile 
generate  high  thrust  force  with  near-zero  heat  dissipation. 

•  The  traveling-wave  constant-flux  concept  provides  the  required  thrust  with  near-zero 
heat  dissipation  in  the  projectile.  This  eliminates  the  inherent  speed  limitation  of 
induction-based  launchers. 

•  Heat  dissipation  in  the  launcher  is  minimal  and  distributed,  resulting  in  higher 
reliability  and  survivability  of  the  launch  tube. 

•  Low-pressure  helium  atmosphere  in  launcher  reduces  aerodynamic  disturbances  and 
aerothermal  heating  during  launch  to  negligible  levels. 

•  Projectile  stabilization  is  aided  by  spinning  in  a  low-pressure  environment  during 
launch.  This  contributes  to  improved  range  and  accuracy. 

•  The  proposed  concept  can  be  used  in  a  wide  range  of  applications  from  kinetic  kill 
weapons  to  direct  LEO  insertion  of  small  payloads. 

•  The  proposed  concept  is  highly  scalable  to  various  payload  sizes  and  launch  speeds. 

•  The  acceleration  profile  can  be  adjusted  according  to  mission  requirements. 

For  these  reasons,  the  proposed  Phase  2  project  constitutes  a  major  step  towards  the 
realization  of  reliable,  responsive,  high-performance  hypersonic  electromagnetic  launch 
systems. 
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1.  Introduction 


Electromagnetic  (EM)  launching  has  been  pursued  in  the  last  three  decades  as  an  efficient 
and  economical  alternative  to  chemical  propulsion  for  military  and  aerospace  applications. 
Although  several  different  concepts  have  been  proposed,  a  number  of  unresolved  technical 
challenges  have  prevented  it  from  becoming  a  reality.  Broadly  speaking,  the  critical  limitations 
of  EM  launching  have  been  (i)  large  heat  dissipation,  leading  to  potential  damage  to  both 
payload  and  launcher  when  reaching  high  speeds,  (ii)  complex  wear  and  contact  phenomena 
resulting  from  mechanical  contact  at  high  speeds,  (iii)  the  lack  of  an  efficient  method  to 
generate  current  in  the  projectile  without  excessive  heating  of  the  payload. 

This  paper  presents  the  constant-flux  synchronous  motor  (CFSM),  a  novel  concept  that 
addresses  most  of  these  critical  problems.  The  CFSM  launcher  is  based  on  a  constant  magnetic 
flux  zone  that  moves  with  the  projectile  in  a  traveling  wave.  Constant-flux  operation  eliminates 
inductive  heating  of  the  projectile.  Superconducting  rings  operating  in  persistent  mode  are  used 
as  current  source  in  the  projectile,  eliminating  the  need  of  feeding  the  launch  vehicle  through 
brushes,  contact  armatures  or  AC  induction.  The  resulting  contact-free,  self  centering 
suspension  and  propulsion  dramatically  improves  survivability  (wear  and  tear)  of  the  launcher 
since  heat  dissipation  in  the  track  during  launch  is  relatively  small  and  distributed  over  the 
length  of  the  launcher.  A  conical/cylindrical  projectile  design  allows  stabilization  by  spinning, 
introduced  at  the  home  position  prior  to  launch-  this  stabilizes  the  projectile  both  during  the 
acceleration  phase  and  during  the  coasting  flight  through  the  atmosphere.  There  is  no  need  for 
detachable  sabot  or  other  parasitic  mass  -  the  effective  launch  mass  can  be  100  %.  For  these 
reasons,  the  CFSM  can  be  a  breakthrough  concept  for  both  direct  EM  launch  to  orbit  and  EM 
kinetic  kill  weapons. 


2.  Development  of  Key  Specifications 

A  detailed  MathCad  program  (see  attachment)  has  been  developed  to  investigate  the 
potential  parameter  space  of  the  proposed  launch  system  and  to  develop  baseline  specifications 
for  the  current  project.  MathCad  was  chosen  since  such  programs  can  be  (almost)  read  like 
normal  documents  and  the  performed  calculations  can  be  easily  followed.  All  quantities  in  the 
MathCad  program  have  physical  dimensions,  which  allows  for  a  simple  cross  check  of  the 
results.  The  complete  program  “Estimate  Electromagnetic  Launch  Assist”  is  given  in  the 
Appendix.  The  current  baseline  parameter  set  is  shown  in  Table  1: 

The  drag  coefficient  of  a  moving  projectile  changes  with  the  Mach  number  showing  a 
significant  enhancement  of  drag  near  Mach  equal  to  1 .  A  typical  profile  as  used  in  the  program 
is  shown  in.  A  more  precise  determination  of  the  aerodynamic  effects  is  given  in  section  4. 

The  electromagnetic  drag  depends  on  the  details  of  the  magnetic  levitation  system,  which  is 
currently  not  determined.  We  have  therefore  conservatively  assumed  that  5%  of  the 
acceleration  force  (without  aerodynamic  drag)  is  needed  to  overcome  this  drag.  The  actual 
value  is  likely  to  be  lower. 
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Table  1:  Input  parameters  for  launch  assist  calculations.  Assumed  base  line  parameters 

The  calculated  performance  values  and  requirements  for  the  chosen  input  parameters  of  Table 
1  are  shown  in 
Table  2. 


Parameter 

Unit 

Value 

Total  mass  to  accelerate 

kg 

105 

T rack  length  for  acceleration 

m 

499.7 

Acceleration  force 

kN 

5410 

Drag  force  in  air 

kN 

1570 

Drag  force  in  helium 

kN 

219 

Total  propulsion  force  in  air 

MN 

6.97 

Total  propulsion  force  in  helium 

MN 

5.63 

Kinetic  energy 

Mjoule 

2570 

Drag  energy  in  air 

Gjoule 

0.39 

Drag  energy  in  helium 

Gjoule 

0.05 

Total  energy  in  air 

Gjoule 

2.96 

Total  energy  in  helium 

Gjoule 

2.63 

Peak  power  in  air 

Gwatt 

48.82 

Peak  power  in  helium 

Gwatt 

39.38 

Table  2:  Calculated  requirements  of  EMLA  system  for  base  line  parameters 
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Drag  Coefficient  versus  velocity 


Figure  1 :  Drag  coefficient  versus  velocity  for  a  drag  coefficient  of  1.  It  is  assumed  for  the  initial 
estimates  that  the  shown  dependence  describes  the  dependence  of  drag  coefficient  on  velocity 
for  our  system. 


2.1  Ambient  Atmosphere 

The  MathCad  program  considers  different  gases  as  the  ambient  atmosphere,  which  are 
described  by  their  density  and  their  specific  heat  ratio  y  =  cp/cv.  For  the  results  shown  in 

Table  2  it  is  assumed  that  the  ambient  atmosphere  consists  of  air  or  helium  at  normal 
pressure  (1  atm).  Different  pressures  can  be  chosen  by  varying  the  density  p  [kg/m3]  of  the  gas. 

Fdrag  —  Vi*  Cdrag  *  p  *  V2  *A 

Equation  1 

Cdrag :  drag  coefficient 
p:  density  of  gas 

v:  velocity  of  object  (sled) 

A:  cross  section 


Since  the  aerodynamic  drag  is  a  significant  contribution  to  the  required  acceleration  force 
(see  below)  it  would  be  of  interest  to  reduce  the  pressure  of  the  atmosphere  surrounding  the 
accelerating  sled.  Also  this  is  not  ruled  out  at  the  present  state  of  the  project,  it  could  be  a 
significant  complication.  The  transition  of  the  sled  form  a  reduced  pressure  to  the  normal 
atmospheric  pressure  at  the  end  of  the  track  could  cause  a  significant  shock  and  has  to  be 
studied  in  more  detail. 

2.2  Track  Length  and  Track  Inclination 

The  track  length  required  for  the  EMLA  system  depends  on  the  launch  speed  and  the 
assumed  acceleration.  The  launch  speed  for  the  proposed  system  is  assumed  to  be  7000  m/sec, 
which  is  about  Mach  20  in  air  and  Mach  7  in  helium.  As  shown  in  Table  1  a  track  length  of 
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500  m  is  assumed.  As  shown  in  our  recent  publication1,  the  track  has  to  be  inclined  at  the  track 
end  to  about  45  degrees  or  more  for  an  efficient  launch.  This  could  be  achieved  by  a  straight 
track  with  constant  inclination  or  a  track  that  starts  horizontally  and  bends  upwards  to  the 
required  angle  at  the  end  of  the  track.  However,  as  shown  below  (Section  3),  the  centrifugal 
forces  acting  on  a  bend  track  would  be  so  large  that  the  solution  would  be  difficult  to 
implement  and  would  compromise  the  overall  system  reliability.  A  track  with  constant 
inclination  is  therefore  chosen  for  the  baseline  design. 

A  lower  limit  of  the  track  length  is  imposed  by  the  increased  acceleration  needed  for 
shorter  tracks.  The  sled  acceleration  will  be  supplied  by  a  linear  drive  system  (TBD),  which 
becomes  increasingly  demanding  for  higher  accelerations.  Also  the  peak  power  requirements 
for  the  linear  drive  system  increases  with  acceleration  and  become  increasingly  difficult  to 
achieve.  As  a  compromise  between  track  length  and  acceleration  we  have  chosen  to  limit  the 
acceleration  to  about  5000g  at  the  end  of  the  track. 

2.3  Deceleration  Section 

A  deceleration  section  would  have  the  advantage  that  a  launch  could  be  interrupted  at  any 
time,  which  would  improve  the  system  reliability.  Such  an  additional  track  section  would  also 
(in  principle)  allow  to  decelerate  the  sled  after  the  projectile  has  been  launched  and  to  reuse  the 
sled  for  another  launch.  However,  the  deceleration  section  would  require  a  significant  increase 
in  track  length  and  cost.  Also  an  additional  linear  drive  system  would  be  needed,  and  the  track 
construction  would  significantly  increase  in  complexity  since  the  launch  point  will  be  at  a 
significant  height  due  to  the  track  inclination. 

With  our  present  understanding  of  the  EMLA  system  we  think  that  the  track  solution 
without  deceleration  section  is  the  preferred  option.  The  cost  of  the  track  with  the  linear  drive 
system  is  significantly  lower  without  the  deceleration,  and  it  is  not  clear  to  which  extend  the 
sled  would  be  reusable  after  a  launch.  The  large  forces  and  surface  temperatures  (see  below) 
experienced  by  the  sled  during  the  acceleration  could  easily  lead  to  some  damage  and  a 
significant  amount  of  inspection  and  repair  work  might  be  required  before  another  launch. 

2.4  Base  Line  Parameters  and  Subsystem  Specifications 

The  subsystem  specifications  developed  with  the  help  of  the  described  MathCad  program 
are  regarded  as  the  design  goals  for  the  various  subsystems  of  the  EMLA  system.  The  detailed 
analysis  performed  in  the  different  project  tasks  might  lead  to  modifications. 

Key  design  requirements  determined  with  the  MathCad  program  are  as  follows: 

The  aerodynamic  drag  acting  on  the  sled  is  significantly  stronger  for  air  than  for  helium. 
The  dependence  for  air  and  helium  at  1  atm  are  shown  in  Figure  2.  The  values  shown  in  the 


1  “Rocket  Performance  Analysis  Using  Electromagnetic  Launch  Assist”,  AIAA  Aerospace 
Sciences  Conference,  Reno  2005,  A.  Uranga,  D.  Kirk  and  H.  Guiterrez,  FIT,R.B.  Meinke, 
AML,K.  Barker,  Schafer  Corporation 
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figure  are  based  on  the  assumed  sled  cross  section  of  0.5  m2  and  a  drag  coefficient  of  0.1  near 
Mach  0. 


velocity  [m/sec] 

Figure  2:  Ratio  of  drag  to  propulsion  force  as  a  function  of  sled  velocity  for  air  (red  curve)  and 
helium  (blue)  curve  based  on  a  track  length  of  500  m  and  a  constant  acceleration  of  5000  g  for  a 
pressure  of  1  atm. 

For  the  assumed  constant  acceleration  of  5000  g  the  force  necessary  to  overcome  the 
aerodynamic  drag  is  about  30%  of  the  force  needed  for  the  acceleration  alone  (without 
aerodynamic  drag)  at  the  end  of  the  track.  The  equivalent  value  for  helium  is  significantly 
lower  and  amounts  to  less  than  4%.  Since  the  acceleration  force  is  constant,  the  curves  also 
reflect  the  increasing  importance  of  the  aerodynamic  drag  with  increasing  speed.  A  detailed 
analysis  of  the  aerodynamic  effects  is  performed  in  Task  3  (see  below). 

The  importance  of  the  aerodynamic  drag  on  the  acceleration  force  has  been  studied  for 
various  acceleration  profiles.  In  air  for  constant  accelerations  larger  than  1400  g  the  plotted 
ratio  from  Figure  2  is  less  than  1,  for  accelerations  below  1400  g  the  aerodynamic  drag 
becomes  more  important  than  the  acceleration  force.  For  helium  this  crossover  happens  for 
acceleration  of  about  200  g. 

It  is  also  of  interest  to  point  out  that  the  peak  in  the  drag  coefficient  near  Mach  1,  as  shown 
in  Figure  1,  has  very  little  effect  on  the  shown  drag  profile.  Only  for  launch  velocities  less  than 
500  m/sec  the  drag  peak  becomes  important. 

2.5  Optimization  of  Acceleration  Profile 

Due  to  the  quadratic  dependence  of  the  aerodynamic  drag  on  velocity  (see  Equation  1 )  the 
drag  is  largest  at  the  end  of  the  track.  It  is  therefore  beneficial  to  use  an  acceleration  profile, 
which  is  not  constant,  but  has  the  largest  acceleration  at  the  start  of  the  track  and  tapers  off 
towards  the  end?  This  would  be  of  particular  interest  for  air,  in  which  the  drag  amounts  to 
about  30%  of  the  total  required  force  at  the  end  of  the  track  and  would  even  be  of  increased 
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importance  for  accelerations  smaller  than  5000  g.  As  part  of  the  described  MathCad  program 
this  optimization  is  performed. 

For  reasons  of  simplicity  a  linear  acceleration  profile  is  assumed,  which  starts  with  an 
initial  acceleration  bstart  and  falls  linearly  with  track  length  during  the  whole  cycle  as  shown  in 
Equation  2. 

Equation  2  acc(s)  =  bstart  +  trio  *  s 

s:  track  length  or  position  on  the  track 

m0:  slope  of  acceleration  profile 

bstart:  acceleration  at  s  =  0 

For  the  real  system  the  acceleration  profile  will  most  likely  contain  a  short  ramp  from  b  =  0 
to  the  initial  value  bstart,  but  this  has  has  no  significant  effect  on  the  performed  optimization. 
Assuming  the  linear  profile  yields  the  following  differential  equation  of  motion: 


Equation  3 


bstart  +  nio-x  =  —  x(t) 
dt 


The  differential,  describing  the  acceleration  profile,  can  be  solved  analytically  for  velocity 
and  position  of  the  sled  as  a  function  of  time  t.  The  following  solutions  are  obtained: 


Equation  4 


Using  Equation  4  and  varying  the  parameters  bstart  and  mo,  the  acceleration  profile  has  been 
optimized  to  make  the  required  propulsion  force  independent  of  track  position.  The  solution  for 
air  is  shown  in  Figure  3  (solid  red  curve). 

The  optimized  acceleration  profile  reduces  the  force  required  from  the  linear  drive  system, 
which  otherwise  rises  towards  the  end  of  the  track.  Also  the  peak  power  requirements  are 
mitigated  by  this  approach.  It  might  even  be  advantageous  to  choose  a  further  reduction  in  the 
acceleration  towards  the  end  of  the  track,  since  it  will  be  more  difficult  to  supply  the  necessary 
acceleration  force  at  the  higher  speeds. 
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Total  Propulsion  forces  in  AIR 


Figure  3:  Propulsion  force  as  a  function  of  track  position  in  air  for  a  constant  acceleration  of 
5000  g  (dashed  curve)  and  a  linear  falling  acceleration  profile  with  bstart  =  5550  g  and  an  end 
acceleration  of  4100  g. 


The  corresponding  propulsion  force  versus  track  position  for  the  same  acceleration  profile 
(bstart  =  5500  g  and  end  acceleration  of  4100  g)  but  for  helium  is  shown  in  Figure  4. 


Figure  4:  Propulsion  force  as  a  function  of  track  position  in  helium  for  a  constant  acceleration 
of  5000  g  (dashed  curve)  and  a  linear  falling  acceleration  profile  with  bstart  =  5550  g  and  an  end 
acceleration  of  4100  g. 

The  reduced  force  at  the  end  of  the  track  (shown  in  Figure  4)  leads  to  a  reduction  in  peak 
power  and  would  facilitate  the  system  design. 


For  Further  details  of  the  MathCad  program  and  the  performed  calculations  see  Appendix  A. 
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3.  Develop  Guideway  Concept 

As  shown  in  the  description  of  Table  1  a  baseline  value  of  500  m  has  been  chosen  for  the 
track  length.  To  avoid  a  complex  control  mechanism  in  the  projectile,  it  has  to  be  released 
under  a  significant  angle  between  25  -  45  degrees  or  higher  (the  optimum  launch  angle  is  still 
TBD).  The  launch  angle  can  be  implemented  by  using  a  sloped  track  or  having  a  modest  bend 
in  the  track  at  the  end  or  over  its  full  length.  A  comparison  of  these  two  cases  is  shown  below. 

A  conceptual  design  for  a  bend  track  is  shown  in  Figure  5.  To  minimize  the  centrifugal 
forces,  the  track  is  curved  over  its  full  length  and  the  circular  arc,  which  starts  horizontally,  is 
truncated  when  it  reaches  45  degrees.  The  arc  length  was  chosen  to  be  the  baseline  track  length 
of  500  m.  The  resulting  radius  is  637  m.  The  track  could  be  fully  underground  to  avoid  any 
above  ground  structures.  The  required  tunnel  would  start  at  a  depth  of  1 87  m. 


tunnel,  such  forces  could  lead  to  movements  of  the  track  during  a  launch.  Under  these 
circumstances  it  might  be  necessary  to  survey  and  realign  the  track  after  each  launch.  For  a 
system  requiring  high  reliability  and  responsiveness  the  curved  solution  does  not  seem 
appropriate  and  a  straight  track  is  preferred.  A  conceptual  solution  is  shown  in  Figure  7. 


Equivalent  Weight  of  Centrifugal  Force 


Figure  6:  Centrifugal  forces  acting  on  the  track  as  a  function  of  sled  velocity.  The  forces  have 
been  converted  to  equivalent  weights  (shown  in  tones)  by  dividing  the  force  values  by  g  =  9.8 
m/sec2. 

To  limit  the  height  of  the  necessary  support  structure  above  ground  the  track  is  inserted 
over  half  of  the  length  into  an  underground  tunnel.  For  the  baseline  track  length  of  500  m  the 
underground  section  would  start  177  m  below  the  surface. 

Tunnels  of  the  required  length  with  the  necessary  underground  facility  do  not  constitute  a 
major  technical  challenge.  Existing  tunnels  for  high  energy  physics  accelerators  (e.g.  the  Large 
Hadron  Collider  (LHC)  tunnel  at  the  European  Center  for  Nuclear  Research,  Geneva, 
Switzerland  is  more  than  30  km  in  length  with  much  larger  underground  facilities). 

A  conceptual  layout  of  the  required  launch  facility  is  shown  in  Figure  8. 

It  is  of  interest  to  point  out  that  the  proposed  launch  concept  would  also  be  applicable  to  a 
canon-like  launch  system  in  which  the  launch  speed  is  reduced  to  about  Mach  2.  Such  a  launch 
velocity  should  be  sufficient  for  shooting  the  sled  with  payload  for  a  few  hundred  miles  (TBD). 
The  centrifugal  forces  for  such  a  system  would  be  rather  modest  (equivalent  weight  of  about  15 
tonnes  for  same  track  length  and  radius),  and  the  bend  track  emerging  from  an  underground 
facility  or  a  even  a  ship  might  be  an  attractive  solution. 
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4.  Description  of  the  CFSM  Launch  Tube  Concept 

The  rail  gun  is  the  EM  launch  technology  that  has  received  the  most  attention  in  recent 
years.  Rail  guns  have  several  attractive  characteristics  but  also  some  unresolved  challenges, 
mostly  given  by  the  complex  contact  phenomena  between  rails  and  armature  at  very  high 
speeds,  as  well  as  the  significant  heat  dissipation  associated  to  each  launch,  which  results  in 
damage  to  the  gun's  insulation  that  substantially  limits  the  rate  of  fire  and  the  overall 
survivability  of  the  gun.  Coaxial  induction  launchers  have  some  desirable  advantages  over 
other  EM  launch  methods  in  that  no  direct  electrical  contact  is  required,  but  pulsed  induction 
motors  or  coil  guns  that  can  provide  the  required  thrust  force  for  hypervelocity  applications 
would  need  to  induce  very  large  currents  in  the  projectile  -  the  resulting  ohmic  losses  would 
lead  to  unmanageable  heat  loads  and  thermal  destruction  of  the  projectile. 

In  order  to  overcome  the  outlined  difficulties  with  existing  propulsion  systems  a  novel 
concept  called  “Constant  Flux  Synchronous  Motor”  (CFSM)  has  been  developed  under  this 
grant.  This  motor  concept  delivers  the  required  propulsion  forces  without  inducing  any 
significant  eddy  currents  in  neither  projectile  nor  launch  vehicle.  To  avoid  inducing  significant 
eddy  currents,  the  projectile  coils  have  to  see  an  almost  constant  magnetic  flux  during  launch,  a 
requirement  that  has  never  been  met  in  any  existing  system.  The  CSFM  launch  system  consists 
of  the  following  components: 

a)  Normal  conducting  coils  on  the  stationary  track,  pulsed  with  30  to  100  kA  during  launch 
depending  on  the  number  of  turns  in  the  windings.  The  pulses  are  synchronized  with  the 
moving  projectile  by  an  appropriate  hardware  trigger.  Near-constant  flux  seen  by  the  projectile 
is  essential  to  the  CFSM  concept.  A  change  in  flux  in  the  projectile's  coils  would  induce  a 
voltage  that  would  create  eddy  currents  in  the  coil's  substrate.  The  resulting  ohmic  heating 
could  quench  or  even  destroy  the  projectile's  coils. 

b)  Superconducting  coils  on  the  projectile,  operated  in  persistent  mode  (30  to  50  kA, 
depending  on  the  number  of  turns).  The  coils  are  inductively  charged  prior  to  launch  -  no 
additional  power  is  needed  to  maintain  the  DC  currents  in  the  projectile  coils. 

c)  Independent  pulsed  power  supplies  connected  to  the  launch  track  coils.  Current  in  the 
launcher  coils  follow  programmed  pulse  shapes  (rise  time,  fall  time  and  pulse  duration) 
triggered  by  the  position  of  the  projectile  along  the  track.  The  shape  of  the  current  pulses  is 
calculated  off-line  -  they  depend  on  the  instantaneous  geometric  relationship  between  track  and 
projectile  coils.  The  pulsed  power  supplies  can  be  implemented  as  2  or  3  combined  Marx 
generators  charged  with  individually  controlled  voltages  and  discharged  with  individually 
controlled  delays  for  each  track  coil,  in  a  way  that  pulse  shapes  can  be  adjusted  to  the  required 
waveforms. 

d)  An  energy  storage  system  with  sufficient  storage  capacity  and  power  rating.  There  are 
several  qualified  technologies  that  could  be  used  (e.g.,  batteries,  capacitors  or  superconducting 
magnetic  energy  storage  devices).  Since  the  heat  dissipation  is  substantially  smaller  than  in 
other  EM  technologies,  and  the  effective  mass  fraction  can  be  100%,  the  energy  and  power 
requirements  are  advantageous  compared  to  other  EM  launch  technologies. 
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Figure  9: 

A  longitudinal  section  of  the  CFSM  launcher.  PB:  Basic  cell  length  of  CFSM  system  =  2  MxPT  , 
for  some  integer  M.  The  distance  between  the  rear  coil  TR1  and  the  corresponding  front  coil  TF1 
is  MxPT.  PP:  distance  between  coils  in  the  projectile,  PT:  distance  between  coils  in  the  track. 
Super-indexes  (p)  and  (n)  denote  coils  of  the  previous  and  next  track  sections  to  be  fired. 


Conventional  coaxial  launchers  rely  on  either  external  power  supplies  that  feed  the 
projectile  through  brushes,  or  induced  eddy  currents.  These  technologies  have  inherent 
limitations  in  the  maximum  speeds  they  can  provide.  In  a  DC  coaxial  launcher  this  is  given  by 
the  wear  and  contact  phenomena  in  the  projectile's  brushes,  as  well  as  ohmic  heating  of  the 
projectile.  In  an  induction  launcher,  inductive  heating  of  the  projectile  is  a  limiting  factor. 

In  a  conventional  pulsed  coil  gun,  the  change  in  speed  of  the  projectile  during  launch  is 
accommodated  by  increasing  either  pole  pitch  or  input  frequency  or  both  along  the  track.  In  the 
CFSM,  the  sequence  of  pulses  supplied  to  the  track  coils  forms  a  traveling  wave  that 
determines  the  speed  of  the  projectile.  Since  there  are  no  iron  cores  in  the  system,  the  pulse 
shapes  required  to  produce  constant  flux  as  seen  by  the  projectile  coils  depend  solely  on  the 
geometric  relationships  of  track  and  projectile  and  can  be  analytically  calculated  off-line. 


4.1  CFSM  Launch  Track  Layout 

The  interaction  between  the  track  and  projectile  coils  provide  the  required  propulsion 
forces  as  well  as  self-centering  suspension  forces.  Projectile  stabilization  during  launch  is 
achieved  both  by  the  self-centering  action  of  the  electromagnetic  forces  and  by  spinning  the 
projectile  prior  to  launch.  The  use  of  a  low  pressure  environment  in  the  track  will  minimize 
drag  and  heat  transfer  during  launch.  The  schematic  layout  of  the  CFSM  launcher  is  shown  in 
Figure  9.  The  basic  cell  of  the  launch  tube  consists  of  M  "rear"  coils  ( Trj ,  Tr2,  ...  Trm)  and  M 
"front"  coils  (7>/,  TF2,  ...  7>w).  The  "rear"  coils  of  such  a  section  are  pulsed  with  opposite 
current  direction  than  the  "front"  coils.  The  distance  between  track  coils  Pt  is  the  distance 

where -  (the  gradient  of  M  with  respect  to  x)  is  maximum.  This  is  calculated  off-line  and 

dx 

depends  only  on  the  launcher  geometry. 
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The  force  between  the  projectile  coil  P  and  the  track  coil  T  is  given  by:  F  -  ipiT - .  The 

dx 

force  is  attractive  when  both  currents  are  in  the  same  direction  and  repelling  when  the  currents 
have  opposite  directions.  While  the  projectile  coils  are  repelled  from  the  nearest  (active)  track 
coil,  the  current  pulse  in  the  track  coil  increases  in  such  a  way  that  it  keeps  the  magnetic  flux 
seen  by  the  projectile  approximately  constant  until  the  projectile  reaches  the  next  track  coil.  At 
this  moment  in  time  the  current  pulse  in  the  active  track  coil  rapidly  goes  to  zero  and  the  next 
track  coil  is  pulsed  in  such  a  way  that  the  flux  through  the  projectile  stays  (approximately) 
constant.  Synchronization  is  achieved  by  hard-wiring  the  firing  circuits  of  individual  energy 
storage  elements  to  sensors  that  detect  the  position  of  the  projectile  along  the  launch  tube. 

In  Figure  9: 

•  The  projectile  coil  currents  ipR  and  ipF  are  constant  (persistent  mode  superconducting 
currents)  and  both  have  the  same  current  direction. 

•  Each  track  coil  is  connected  to  an  independent  pulse  power  supply  which  is  able  to  supply 
a  pre-determine  pulse  shape.  Only  two  track  coils  are  ON  at  any  given  time.  In  the  picture,  Tri 
and  Tfi  (in  red)  are  currently  being  pulsed.  Triggering  the  track  coils  is  hardware-controlled  by 
proximity  switches  that  detect  the  presence  of  the  projectile. 

•  The  "rear"  pulse  im(t)  starts  when  the  projectile  coil  Pr  reaches  x  -  -PT  and  ends  when  Pr 
reaches  x  =  0.  The  coordinate  x  is  relative  to  the  coordinate  system  that  moves  with  the 
projectile.  Similarly,  the  "front"  pulse  iFi(t)  is  only  ON  while  the  front  projectile  coil  PF  is 
traveling  through  the  front  shaded  area.  The  current  pulses  iR!(t)  and  iFi(t)  are  calculated  off¬ 
line  as  described  in  the  next  Section.  The  rear  coil  pulls  the  projectile  to  the  right  (currents  in 
the  same  direction)  while  the  front  coil  pushes  it  to  the  right  (currents  have  opposite 
directions). 

•  After  a  section  is  fired,  the  next  section  starts  ( Tr /n),  Tr2("\  ...).  Notice  that  the  "front"  coils 
of  the  previous  section  fired  are  placed  in  close  proximity  to  the  "rear"  coils  of  the  current 
section  (Tri  with  Tfip\  Tr2  with  Tf2(p),  and  also  Tr/"*, with  Tfi,  TR2(n)with  Tf2,  etc.  ). 

•  For  a  slender  projectile  (i.e.  length  »  diameter),  it  is  assumed  that  the  effect  of  a  "front" 
track  pulse  on  the  rear  shaded  section  is  negligible.  Similarly,  the  effect  of  the  "rear"  track 
pulse  on  the  front  shaded  section  is  negligible  too. 

•  The  mutual  inductance  M(x)  between  the  projectile  coil  and  the  track  coil  that  are 
interacting  at  any  given  moment  (in  the  picture,  PR  with  Tri,  and  PF  with  Tfi)  is  a  function  of 
the  instantaneous  axial  position  V  and  the  basic  launcher  geometry  (launcher  diameter, 
projectile  diameter,  and  air  gap).  Since  there  are  no  iron  cores,  M(x)  is  easily  and  accurately 
calculated  off-line. 

Coils  Tri  and  Tfi  are  fired  at  the  same  time,  when  x  =  -Pr-  When  x  =  0,  Tri  and  Tfi  are 
turned  OFF.  The  origin  of  the  coordinate  system  "X"  moves  to  Tri.  At  the  same  time,  the 
pulses  in  Tr2  and  Tf2  start.  Again,  the  "rear"  currents  are  fired  in  the  same  direction  as  the 
projectile  current,  whereas  the  "front"  coils  are  fired  in  opposite  direction.  This  way,  the  rear 
coil  is  always  pulling  the  projectile  forward  while  the  front  coil  is  pushing.  Tr2  and  Tf2  are 
turned  OFF  when  the  local  coordinate  x  reaches  zero.  The  sequence  is  repeated  until  the  last 
coils  of  the  section  Trm  and  Tfm  are  fired.  The  next  section  (Tr/”*  and  Tf/”\  Tr and  7>/^, 
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...)  fires  in  similar  manner.  Notice  that  the  coils  in  the  "front"  section  always  fire  in  direction 
opposite  to  the  projectile  coil  currents,  and  so  on.  If  a  given  pulse  shape  is  difficult  to  achieve 
with  a  single  pulse  source  or  Marx  generator,  two  or  three  can  be  operated  in  parallel  to  obtain 
the  required  pulse  shape  by  superposition.  Some  pulse  characteristics  (i.e.  rise  time)  are  given 
by  track  parameters  such  as  number  of  turns  (i.e.,  inductance)  of  a  given  coil. 

4.2  Analytical  Overview  of  the  CFSM  Concept 

As  described  in  the  previous  section,  the  only  track  coils  active  at  any  given  time  are  the 
two  located  immediately  ahead  (Tri)  and  immediately  behind  (7>y)  of  the  projectile  rear  and 
front  coils,  respectively.  The  corresponding  current  pulses  last  only  while  coils  Pr  and  Pf  are 
moving  within  their  corresponding  shaded  areas.  The  time  taken  to  travel  the  shaded 
subsection  Pt  is  tpj. 

The  magnetic  field  and  flux  acting  on  the  projectile  coils  Pr  and  Pf  while  traveling  within 
the  shaded  areas  (see  Figure  9)  will  be  considered.  A  simpler  case  will  be  analyzed  first:  pitch 
and  yaw  angles  are  assumed  to  be  small,  and  the  centers  of  coils  PR  and  Pf  are  considered  near 
the  axis  of  symmetry  of  the  launch  tube.  In  this  case,  the  magnetic  field  B  produced  by  the 
track  coils  Tri  and  Tpi  on  any  circular  concentric  path  within  the  surfaces  enclosed  by  Pr  and 
Pf  has  cylindrical  symmetry  (i.e.,  B  is  the  same  for  all  points  along  any  concentric  circular  path 
within  the  surfaces  enclosed  by  Pr  and  Pf).  The  flux  ^caused  by  the  track  coil  Tri  on  the 
surface  enclosed  by  PR  can  therefore  be  calculated  using  circular  differential  elements  as 
shown  in  Figure  10. 


Pr  TRi 


Figure  10:  Magnetic  Flux  through  a  Projectile  Coil  section. 

The  flux  ^caused  by  the  track  coil  Tri  on  the  surface  enclosed  by  PR  is: 

N 

tj>  =  \BdA  =  ^B(27±Ar)Ar ,  where  B  is  some  nonlinear  function  of  the  geometry  and  the 

A  k= 0 

current  />/,  B  =  iRlg(r,x) ,  and  r  =  kAr  is  the  discretized  version  of  the  continuous  variable  r. 
The  index  k  numbers  the  concentric  rings  used  in  this  numerical  approximation.  The  field  of 
the  projectile  coil  Pr  is  not  included  on  the  flux  calculation  since  it’s  a  stationary  constant  field 

that  travels  with  the  projectile  and  therefore  causes  no  induction  e  =  on  the  projectile’s 

coil.  The  condition  for  near-constant  flux  on  the  projectile  coil  Pr  traveling  in  the  rear  shaded 
area  of  Figure  1  can  now  be  written  as: 
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(1) 


N  N 

2/r^  iRi  (t)g(kAr,  x)(kAr)Ar  +  iF]  (t)g(kAr,  x  +  MPT  )(kAr)Ar  =  O  0 
*= o  *=o 

which  includes  the  effect  of  the  two  track  coils  being  pulsed  (7>/  and  7>/,  with  currents  iRi(t) 
and  ipi(tj)  on  the  rear  projectile  coil  Pr.  M*Pt  is  the  distance  between  coils  Tri  and  7>;  (the 
length  of  a  basic  track  section  is  2M*PT)  and  O0is  some  characteristic  flux  to  be  kept  near 
constant.  For  a  slender  projectile  (i.e.,  the  projectile  radius  Rp  is  much  smaller  than  the  length 
of  the  half-section  M*Pt),  this  can  be  simplified  by  neglecting  the  effect  of  the  front  track  coil 
Tfi  (which  is  rather  far  away)  on  the  rear  projectile  coil  Pr. 

2/r^z*!  (t)g(kAr,x)(kAr)Ar  =  O0>  -PT<x<0  (2) 


and  A ,r  =  — ,  where  N  is  an  arbitrarily  large  integer.  This  can  be  regrouped  as: 
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The  time  taken  to  travel  the  shaded  subsection  ( tpp )  is  discretized  as  t  =  jAt ,  where 

At  =  —  and  M,  is  some  (large)  integer ,j=0  ...  M, .  The  required  current  pulse  iRI(t)  at  time 
M, 


t  =  jAt  ,  must  meet  the  near-constant  flux  condition  (3): 
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where  x  =  x(  jAt)  is  the  position  of  the  projectile  coil  Pr  relative  to  the  coordinate  system  Tr  at 
the  present  moment.  In  this  expression,  the  required  current  pulse  iRi(t)  does  not  depend  on  the 
integration  variable  k,  and  can  therefore  be  factored  out  of  the  summation  above.  This  yields: 
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for  -PT  <x<  0,  which  illustrates  how  the  desired  current  pulse  hi(t)  is  a  function  of  the 
instantaneous  position  x  =  x(jAt)  ,  and  g(r,x),  that  depends  on  the  geometric  relationship  of 
track  and  coil.  The  current  iR/(t)  and  the  position  x  =  x(jAt)  are  related  in  this  model  through 
Newton’s  second  law  and  the  electromagnetic  thrust  forces: 
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for  -PT  <x<0,  where  m  is  the  projectile’s  mass,  ipn  is  the  projectile  coils  current  (a  constant) 

and  the  gradient  - =  / (x)  can  be  calculated  off-line  and  interpolated  as  a  function  of  the 

dx 

projectile’s  position  x.  In  the  launch  tube’s  low-pressure  atmosphere  the  aerodynamic  drag  can 
be  neglected  and  therefore  the  total  force  in  the  x  direction  is  determined  by  the 
electromagnetic  thrust  forces.  The  dynamic  dependence  of  current,  force  and  position  can  be 
described  by  converting  the  algebraic  condition  (5)  into  a  differential  condition  in  time.  Taking 
the  derivative  of  Equation  (3)  with  respect  to  time,  yields,  for  -PT  <  x  <  0 : 
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where  r  =  k A  similar  equation  can  be  written  for  the  front  coil  current  iFi(t).  Since  both 

iRi(t)  and  its  time  derivative  are  independent  of  the  summation  variable  k,  this  equation  can  be 
rearranged  as  follows: 
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where  the  projectile  speed  is  large  (x»  r)  and  therefore  the  term  —  r  in  (7)  can  be 

dr 

neglected.  The  corresponding  expression  for  the  front  coil  is: 
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Equations  (6),  (8)  and  (9)  define  the  dynamic  relationship  between  the  track  pulses  ipi(t), 
ifi(t)  and  the  projectile’s  motion.  The  coupled  equations  can  now  be  described  as  a  system  of 
four  nonlinear  ordinary  differential  equations  of  the  four  state  variables: 
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The  coupled  system  of  ordinary  differential  equations  (11)  can  be  readily  solved  using 
standard  numerical  techniques  such  as  the  Runge-Kutta-Fahlberg  or  Adams-Moulton 
algorithms,  starting  at  the  initial  conditions: 

z.  00  =  z2  00=G/ma*»  z3{e)*-PT,  zA(e)*x 0  (12) 


where  £•  is  the  (hardware  given)  shortest  time  in  which  the  energy  storage  devices  can 
discharge  the  maximum  current  peak  imax,  and  x0  is  the  projectile  velocity  at  the  end  of  the 

previous  track  segment  (s  is  typically  very  small  compared  to  the  total  duration  of  the  pulse 
tpr).  The  numerical  solution  of  (11)  subject  to  (12)  ends  when  z3  =  0 ,  yielding  the  current 
waveforms  im(t),  iFl (t)  that  generate  maximum  thrust  in  the  launch  vehicle  while  keeping  a 
near-constant  flux  region  that  travels  synchronously  with  it.  Notice  that  for  the  rear  coil 
solution,  the  initial  condition  is  ipi(s)=  imax  ,  whereas  for  the  front  coil  Ifi(£)~  G  imax  where  0 
<  G  <  1  is  a  constant  to  be  optimized  by  successive  solutions  of  (1 1)  subject  to  (12). 

3.2.1  Computation  of  Mutual  Inductance 

To  estimate  the  electromagnetic  forces  (thrust,  drag  and  guidance)  acting  on  the  projectile 
during  launch,  the  conductors  in  the  track  and  projectile  coils  are  first  approximated  by 
infinitely  thin,  circular  filaments.  This  approximation  is  sufficient  at  the  conceptual  design 
stage  and  will  be  improved  during  Phase  2  to  include  the  effect  of  conductor  thickness.  The 
magnetic  field  generated  by  a  circular  filament  conductor  is  given  by: 
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where  and  Bz  are  the  radial  and  axial  component  of  the  magnetic  field,  a  is  the  radius  of  the 
wire  loop,  I  is  the  current  in  the  loop  and  p0  is  the  permeability  of  free  space  (4nx  1 CT7  volt 
sec/A.m).  K  and  E  are  the  first  and  second  kind  of  elliptical  integrals: 
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where  r  and  z  are  the  polar  coordinates  of  any  point  in  space  where  the  magnetic  field  is 
calculated.  Computation  of  the  electromagnetic  forces  requires  an  estimate  of  the  mutual 
inductance  between  the  coils  in  the  track  and  the  corresponding  coils  in  the  projectile,  as 
shown  in  Eq.  (6).  The  basic  computation  of  mutual  inductance  is  for  two  filament  loops  in 
space,  i.e.,  it  is  assumed  that  the  depth  of  penetration  of  the  magnetic  field  in  the  filament  is 
equal  or  greater  than  the  diameter  of  the  filament,  or,  in  other  words,  that  the  current  density  in 
the  filament  is  constant. 


Figure  11:  Mutual  inductance  versus  radial  displacement  for  various  axial  positions  (Z) 
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X-potition  of  the  center  of  the  secondary  coil  (meters) 


Figure  12:  Mutual  inductance  (M)  versus  radial  displacement  for  various  angles  of  rotation  of 

the  secondary 

Two  angles  (Alpha  and  Beta)  define  the  tilt  of  the  moving  coil  relative  to  the  fixed  coil. 
The  coil's  mesh  can  now  be  expressed  relative  to  the  fixed  coil  for  any  rotation  of  the 
secondary.  Rotation  around  the  Z-axis  is  not  considered  due  to  axial  symmetry. 

The  magnetic  field  is  calculated  on  every  point  on  the  secondary's  mesh:  both  radial  (Br) 
and  axial  (Bz)  components  are  found  in  the  fixed  coordinate  system  attached  to  the  primary.  To 
calculate  the  flux,  the  component  of  the  field  perpendicular  to  secondary  surface  is  calculated, 
using  the  inverse  of  the  rotational  matrices  M0i  and  M12.  The  flux  through  the  projectile  coil  is 
defined  as  the  sum  of  all  normal  components  of  the  magnetic  field  expressed  in  the  coordinate 
system  2.  The  flux,  divided  by  the  current  in  the  primary  coil,  gives  the  mutual  inductance. 
Several  geometric  configurations  were  simulated,  using  ^(primary)  =  0.5  m,  Rs  (secondary)  = 
0.2  m.  and  theta  varies  from  0  to  2  n  with  a  step  of  0. 1  rad.  The  radial  coordinate  r  varies  from 
0.01  to  Rs  with  a  step  of  0.01  m.  Figure  11  shows  the  mutual  inductance  as  a  function  of  the 
transverse  motion  (secondary  and  primary  move  in  parallel  planes  for  different  axial  positions). 
Alpha,  Beta  and  the  y  were  set  to  0.  As  the  axial  distance  increases,  the  mutual  inductance 
decreases  as  expected.  For  small  axial  distance  between  primary  and  secondary,  the  mutual 
inductance  has  two  peaks  and  a  local  minimum.  This  behavior  is  related  to  the  local  curvature 
of  the  magnetic  field,  which  is  more  important  at  smaller  axial  distances.  Figure  12  shows  the 
effect  of  the  rotation  angles.  The  axial  position  is  set  to  0.3  m  with  respect  to  the  center  of  the 
primary  coil,  and  a  transverse  motion  in  the  radial  direction  is  simulated.  The  resulting  mutual 
inductance  is  plotted  for  various  rotation  angles. 

As  expected,  the  mutual  inductance  is  higher  when  Alpha  and  Beta  are  equal  to  0  (the  flux 
linkage  is  maximal).  For  Alpha  =  0  and  Beta  =  0,  M  is  near  zero,  except  when  the  secondary  is 
near  the  edges  of  the  primary,  in  which  case  the  magnetic  field  has  a  strong  local  curvature 
around  the  edges,  i.e.,  a  substantial  component  that  generates  flux.  Similarly,  when  Alpha  =  0 
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and  Beta=  n/2.  Depending  on  the  axial  position,  the  mutual  inductance  might  be  negative,  since 
the  flux  depends  on  the  local  curvature  of  the  magnetic  field. 


Figure  13  shows  the  evolution  of  the  mutual  inductance  as  a  function  of  axial  displacement 
between  the  two  coils  for  various  radial  displacements.  For  the  case  shown,  Alpha  and  Beta  are 
zero.  As  expected,  the  mutual  inductance  is  symmetrical  with  respect  to  Z=0  and  maximal 
when  the  two  coils  are  in  the  same  plane.  The  mutual  inductance  increases  with  radial  off-set. 
No  sign  change  of  the  mutual  inductance  is  expected  as  function  of  axial  position. 


Figure  13:  Mutual  inductance  (M)  versus  axial  displacement  (Z),  for  various  center  offsets  d 

between  the  two  loops. 


23 


5.  Persistent  Current  Superconducting  Coils 


As  described  in  Section  4,  the  propulsion  system  requires  secondary  coils  on  the  projectile 
that  interact  with  the  pulsed  coils  on  the  launch  tube  track.  The  operational  current  of  the 
projectile  coils  is  estimated  to  be  30  to  50  kA.  The  mutual  inductance  between  the  track  and 
projectile  coils  has  to  be  as  large  as  possible  to  achieve  the  necessary  propulsion  forces  of 
about  5  MN  (Table  1).  The  required  currents  of  the  projectile  coils  are  sustained  by 
superconducting  ring  coils  mounted  coaxially  on  the  projectile  as  shown  in  Figure  14. 

Superconductivity  is  the  breakthrough  technology  that  allows  large  current  in  the  projectile 
without  brushes,  external  connections  or  the  heat  load  typically  associated  to  large  currents. 
Superconducting  coils  have  practically  zero  resistance  (~  10'9  ohm),  which  is  essential  to  avoid 
excessive  heat  loads  at  the  required  operational  currents  (up  to  50  kA).  These  superconducting 
coils  are  operated  in  a  persistent  mode  and  no  power  supply  is  required  during  launch.  With  an 
estimated  resistance  of  per  coil  of  10‘9  ohm,  the  heat  load  at  50  kA  is  only  2.5  W. 

The  projectile  could  be  launched  as  a  single  object  (no  separation  required  after  launch)  or 
could  detach  after  launch,  if  used  to  deliver  a  payload  such  as  a  laser  propulsion  engine  (see 
Figure  14,  left).  This  flexibility  permits  to  reach  100%  effective  payload  mass  as  desired  in 
kinetic  kill  applications. 

The  superconducting  rings  are  charged  inductively  to  the  required  operational  current  prior 
to  launch.  The  individual  superconducting  wires  or  tapes  are  inserted  in  a  sturdy  metallic 
conduit  filled  with  cryogen  (most  likely  helium  or  hydrogen)  to  cool  the  superconductor  to  its 
required  operational  temperature.  Prior  to  a  launch,  the  coil  conduits  would  be  connected  via 
umbilical  cords  (small  cryogenic  transfer  lines)  to  a  refrigerator,  cooled  to  the  operational 
temperature  and  filled  with  a  certain  amount  of  cryogen.  The  umbilical  cords  are  disconnected 
shortly  before  the  launch,  and  the  superconductor  is  kept  at  the  operational  temperature  by  the 
amount  of  cryogen  inside  of  the  conductor  conduit.  Since  the  launch  duration  is  only  about  0.5 
seconds  (see  Table  1),  aerothermal  heating  is  estimated  to  be  about  50  K  in  the  surface  of  the 
projectile  (see  section  6)  for  a  Mach  20  launch  speed,  with  insignificant  induction  heating 
during  launch  (see  Section  4).  Therefore,  a  small  amount  of  cryogen  is  sufficient  to  keep  the 
conductor  in  the  superconducting  state. 

Under  these  operational  conditions,  the  cryostat  surrounding  these  coils  can  be  simple  and 
inexpensive.  While  cryostats  for  most  applications  require  a  careful  reduction  of  external  heat 
loads,  the  cooling  requirements  for  this  application  are  rather  modest.  The  launch  tube  low- 
pressure  atmosphere  will  function  as  an  insulation  vacuum,  and  very  little  convective  heating 
from  the  surroundings  is  expected. 

Placing  LTS  superconductors  in  metallic  conduits  is  well  established  technology,  widely 
used  in  pulsed  fusion  magnets.  The  direct  contact  between  the  superconductor  and  the  cryogen 
in  this  cable-in-conduit  (CICC)  technology,  enables  that  energy  dissipated  from  pulsed 
magnetic  fields  or  conductor  motion  is  absorbed  by  the  cryogen.  A  typical  superconductor  with 
a  current  carrying  capacity  of  about  50  kA  is  shown  in  Figure  15. 


24 


Figure  14:  left)  Layout  of  projectile  inside  a  section  of  the  launch  tube,  when  used  as  first  stage 
to  launch  a  laser  propulsion  engine,  right)  Projectile  with  superconducting  rings 


Figure  15:  Cable-in  conduit  used  on  the  ITER  CS  Model  Coil.  The  conduit  dimensions  are 
-  50  x  50  mm.  This  CICC  conductor  operated  successfully  at  46  kA  and  13  Tesla,  at  a  magnetic 
field  pulse  rate  of  2  Tesla/sec. 

While  the  specific  heat  of  almost  all  materials  is  extremely  low  at  temperatures  of  about  4.5 
K  (<  1  Joule  /  (kg  K) ),  the  specific  heat  of  helium  peaks  at  this  temperature  and  reaches  values 
of  more  than  8000  Joule  /  (kg  K).  Due  to  its  large  specific  heat  at  these  temperatures,  helium 
can  absorb  significantly  larger  amounts  of  heat  without  significant  temperature  rise. 

The  specific  heat  of  low  temperature  superconductors  at  typical  operational  temperatures  of 
4.5  K  is  also  very  low  -  these  metallic  alloys  cannot  absorb  significant  amounts  of  heat  without 
heating  up  above  their  critical  temperature  and  becoming  normal  conducting.  On  the  other 
hand,  BSCCO,  YBCO  and  MgB2  high  temperature  superconductors  (HTS)  are  already 
available  in  wire  and  tape  form.  If  operated  at  temperatures  below  20  K,  these  conductors  can 
carry  large  currents  at  rather  high  magnetic  fields.  BSCCO-2212,  and  MgB2  operated  at  4.5  K, 
have  a  much  higher  critical  field  than  any  low  temperature  superconductor  (LTS).  While  not 
absolutely  necessary  for  the  proposed  launcher  application,  HTS  conductors  would  offer  a 
much  higher  energy  margin  and  therefore  quench  safety  than  LTS  conductors.  The  amount  of 
heat  that  can  be  accommodated  in  HTS  conductors  can  exceed  those  of  LTS  by  a  factor  of 
several  thousand.  This  is  due  to  the  much  higher  enthalpy  of  materials  at  higher  temperature 
and  the  much  larger  gap  between  operational  temperature  and  critical  temperature.  HTS 
superconductors,  operating  at  20  to  30  K,  therefore  offer  a  much  higher  safety  margin  between 


25 


nominal  operation  and  quench  compared  to  LTS  superconductors  operating  at  liquid  helium 
temperatures  (typically  4.5  K). 

On  the  other  hand,  the  critical  currents  of  HTS  superconductors  (YBCO  and  BSCCO)  can 
be  pressure  sensitive  and  their  operation  under  large  stress  levels  has  not  been  well  established. 
Specific  tests  will  be  required  to  qualify  these  conductors  for  the  launcher  application. 
Appropriate  tests  will  be  proposed  for  a  second  Phase  of  this  project. 

6.  Aerodynamic  and  Aerothermal  Effects 

During  the  Phase  I  effort,  significant  progress  has  been  made  on  understanding  and 
modeling  the  aerodynamics  and  heat  transfer  of  the  sled/projectile  system  while  contained 
within  a  launch  tube.  The  calculations  have  been  performed  using  both  closed-form  analytical 
solutions  as  well  as  using  numerical  calculations. 

Employing  an  evacuated  launch  tube  (~  1/1 00th  of  an  atmosphere)  the  aerodynamic  forces, 
i.e.  lift,  drag  and  moments,  are  within  manageable  levels,  which  can  easily  be  compensated  for 
by  the  electromagnetic  forces.  For  example,  the  lift  and  drag  are  on  the  order  of  1-15  kN  and 
the  moments  at  maximum  misalignment  of  the  projectile  allowed  by  the  launch  tube  are  around 
20  kN  m.  A  parametric  study  has  been  performed  over  a  range  of  Mach  numbers  and  projectile 
orientations.  Since  the  geometry  of  the  sled/projectile  has  changed  to  a  conical/cylindrical 
shape  during  the  course  of  this  investigation,  a  new  analytical  and  computational  study  was 
completed  for  this  shape,  and  the  resultant  lift,  drag  and  moments  being  of  the  same  order  of 
magnitude  as  the  flat  plate  approximation  or  less  due  to  a  three-dimensional  relieving  effect 
associated  with  the  conical  shape.  A  table  of  drag  and  lift  coefficients  is  presented  later  in  this 
section. 

Although  aerodynamic  forces  and  moments  are  easily  overcome  by  magnetic  forces  within 
the  launch  tube  at  relatively  modest  vacuum  levels,  it  has  been  found  that  1/1 00th  of  an 
atmosphere  vacuum  is  essential  to  alleviate  heat  transfer  concerns  within  the  launch  tube. 
Using  stainless  steel  as  the  sled  material,  a  series  of  simulations  were  conducted  to  ensure  that 
the  material  would  not  approach  melting  temperatures.  The  cases  were  performed  by 
evaluating  the  heat  flux  along  the  surface  of  the  projectile  during  a  0.5-second  simulated 
acceleration  from  zero  to  7  km/s.  Even  using  worst  case,  constant  heat  flux  it  was  found  that 
the  surface  temperature  of  the  stainless  steel  body  only  increases  by  about  50  K,  if  the  pressure 
within  the  launch  tube  is  kept  to  1/1 00th  of  an  atmosphere.  To  summarize,  the  preliminary 
studies  have  indicated  that  the  level  of  vacuum  required  is  governed  by  heat  transfer 
considerations  rather  than  the  need  to  minimize  aerodynamics  disturbance  forces  within  the 
launch  tube.  Heat  transfer  within  the  launch  tube  was  also  investigated  using  computational 
methods  which  are  coupled  to  the  flow  solution. 

In  order  to  understand  the  aerodynamics  of  a  projectile  traveling  at  supersonic  and 
hypersonic  speeds  in  the  channel  of  an  electromagnetic  launcher,  numerical  simulations  were 
performed  using  the  Computational  Fluid  Dynamics  (CFD)  commercial  state-of-the-art  code 
FLUENT  version  6.2.  In  this  section,  the  main  results  of  two-dimensional  simulations  are 
presented  and  analyzed. 
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6.1  Analysis  of  a  Simpler  Case:  Diamond-Shaped  Projectile  in  a  Large  Channel 

The  purpose  of  this  section  is  to  provide  a  qualitative  understanding  of  the  key  phenomena 
involved  in  supersonic  flows,  namely  shocks  and  expansion  fans.  Consider  a  projectile  with  a 
diamond  shape,  -0.63  m  long  and  0.04  m  thick,  its  length  and  thickness  being  the  same  as  the 
baseline  sled.  The  walls  of  the  channel  are  placed  sufficiently  far  from  the  projectile  so  that  the 
reflected  shocks  and  expansion  waves  do  not  touch  the  projectile.  The  projectile  is  traveling  at 
Mach  2  in  air  at  1  atm  (101325  Pa)  and  300  K,  and  the  fluid  is  assumed  inviscid  for  this 
preliminary  case.  The  grid  for  the  simulation  is  obtained  through  adaptive  refinement  and  can 
be  seen  in  Figure  16. 

The  velocity,  pressure,  density,  and  temperature  fields  around  the  projectile  are  shown  in 
Figure  17.  An  oblique  shock  appears  at  the  leading  edge  of  the  projectile  across  which  velocity 
decreases  while  pressure,  density,  and  temperature  increase;  a  shock  represents  a  discontinuity 
in  the  flow  properties.  From  the  upper  and  lower  comers  of  the  shape,  an  expansion  fan 
emanates  to  turn  the  flow.  The  flow  properties  change  continuously  across  the  fan:  velocity 
increases  while  pressure,  density,  and  temperature  decrease.  The  drag  force  on  the  projectile  is 
the  result  of  the  difference  in  pressure  between  the  front  half  and  back  half  of  the  shape. 
Because  of  the  expansion  wave,  the  pressure  in  the  front  will  be  higher  than  the  pressure  in  the 
back,  and  hence  a  net  aerodynamic  force  acts  backwards. 

Finally,  as  can  be  seen  in  Figure  17,  both  shocks  and  expansion  waves  are  reflected  by  the 
channel  walls.  The  main  effect  of  the  interaction  between  the  expansion  fan  and  the  shocks  is 
the  bending  of  the  shock  away  from  the  fan.  As  will  be  emphasized  in  the  following  sections, 
even  though  the  bending  angle  is  relatively  small,  this  interaction  can  have  important 
consequences  on  the  overall  aerodynamic  force  on  the  projectile  by  changing  the  location 
where  the  reflected  shock  will  encounter  the  projectile  again. 


Figure  16:  Grid  used  for  the  diamond-shaped  projectile  in  inviscid  flow 
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(a)  Velocity 


(b)  Gauge  pressure 


(c)  Density 


(d)  Temperature 


Figure  17:  Velocity  (a),  pressure  (b),  density  (c),  and  temperature  (d)  fields  for  the  inviscid  flow 
around  a  diamond-shaped  projectile  traveling  in  air  at  Mach  2. 

6.2  Projectile  Traveling  in  a  Narrow  Channel 

The  viscous  supersonic/hypersonic  flow  through  a  very  narrow  gap  between  two  walls  is 
extremely  difficult  to  simulate  mainly  because 

-  the  flow  is  highly  energetic  and  highly  compressible; 

-  each  time  a  shock  reflects  off  a  wall  it  interacts  with  the  boundary  layer,  resulting  in  a 
very  complex  physical  phenomena; 

-  the  turbulence  models  that  need  to  be  used  to  simulate  turbulent  flows  are  accurate  near 
wall  because  they  rely  on  “classical”  boundary  layer  theory,  but  they  are  often 
inadequate  for  use  in  shock-boundary  layer  interactions. 

Because  of  these  issues,  the  numerical  simulation  has  problems  capturing  the  physics  of  the 
flow  and  is  hence  very  unstable  (i.e.  it  tends  to  diverge).  In  the  design  baseline  the  gap  between 
the  sled  and  the  track  is  1cm,  the  sled  being  4  cm  high.  The  supersonic/hypersonic  flow  in  this 
tiny  channel  was  not  possible  to  simulate  accurately.  Therefore,  for  the  computational  fluid 
dynamics  study,  the  gap  was  increased  to  8cm,  which  corresponds  to  a  launcher  channel  of 
10cm  radius.  The  fluid  is  taken  to  be  helium  and  assumed  to  behave  like  an  ideal  gas. 


The  purpose  of  this  analysis  is  to  aid  in  the  development  of  a  simplified  semi-empirical 
model  that  will  describe  the  main  features  of  the  flow  and  predict  the  aerodynamic  drag,  lift, 
and  moment  on  the  projectile  at  any  instant  during  its  flight  inside  the  launcher.  The  starting 
point  of  such  a  model  was  the  analytical  computation  of  the  inviscid  flow  around  the  projectile, 
which  includes  shocks  and  expansion  waves  but  no  reflection  at  the  walls.  From  the  numerical 
simulations  of  both  inviscid  flow  and  viscous  flow  inside  the  channel,  the  importance  of  the 
physical  phenomena  not  captured  in  the  simplified  analytical  model  is  determined.  Corrective 
factors  are  then  introduced  into  the  analytical  model  to  account  for: 


-  reflections  on  the  wall; 

-  interactions  between  shocks  and  expansion  waves; 

-  viscous  effect,  such  as  development  of  a  boundary  layer  on  the  projectile  and  sock¬ 
boundary  layer  interactions. 

These  correction  factors  will  be  applied  to  the  drag,  lift,  and  moment  coefficients,  and  will  be 
functions  of  the  Mach  number,  gap  size  between  the  projectile  and  the  channel  walls,  and 
projectile  angle  of  attack. 

Figure  18  through  Figure  24  show  the  velocity,  pressure,  and  temperature  fields  inside  the 
channel.  Inviscid  cases  are  presented  in  Figure  18  through  Figure  20  and  Figure  25,  while 
viscous  results  can  be  seen  in 
Figure  21  through  Figure  24. 
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Figure  18:  Projectile  traveling  at  Mach  2  through  inviscid  helium  at  1  atm. 
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Figure  19:  Projectile  traveling  at  Mach  5  through  inviscid  helium  at  latm. 
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Figure  20:  Projectile  traveling  at  Mach  7  through  inviscid  helium  at  latm. 
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Figure  21:  Projectile  traveling  at  Mach  2  through  viscous  helium  at  latm. 
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(a)  Closer  view  of  the  leading  edge  (b)  Bloser  view  of  trailing  edge 


(c)  Closer  view  of  the  shock-boundary  layer  interaction 

Figure  22:  Closer  view  on  the  velocity  field  around  a  projectile  traveling  at  Mach  7  through 

inviscid  helium  at  latm. 
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(a)  Closer  view  of  the  leading  edge 


(b)  Closer  view  of  trailing  edge 


(c)  Closer  view  of  the  shock-boundary  layer  interaction 

Figure  23:  Closer  view  on  the  temperature  field  around  a  projectile  traveling  at  Mach  7  through 

inviscid  helium  at  latm. 
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Figure  24:  Projectile  traveling  at  Mach  7  through  inviscid  helium  at  0.001  atm. 
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As  can  be  seen  in  Figure  23  (a),  the  boundary  layer  on  the  projectile  is  very  thin  at  such 
high  speeds  because  the  flow  carries  a  lot  of  inertia.  Figure  23  (c)  shows  the  thin  wake  being 
shed  off  the  back  of  the  projectile.  When  a  shock  is  reflected  by  the  projectile  wall,  it  interacts 
with  the  boundary  layer  and  causes  boundary  layer  separation  as  illustrated  in  Figure  23  (b). 

In  the  viscous  flow  at  Mach  2  of  Figure  24  (a),  one  can  see  a  significant  increase  in 
temperature  close  to  the  walls  caused  by  the  fluid  viscosity  which  dissipates  the  fluid  energy  as 
heat.  Figure  24  (b)  shows  the  high-temperature  flow  being  carried  away  by  the  wake. 
Moreover,  the  highest  temperatures  on  the  sled  are  linked  to  viscous  heating  at  a  shock¬ 
boundary  layer  interaction,  as  can  be  seen  in  Figure  24  (c). 

The  comparison  of  the  inviscid  and  viscous  cases  simulated  shows  that  the  viscous  effects 
increase  the  drag  coefficient  of  the  projectile  by  between  10%  and  15%.  For  instance  the  drag 
coefficient  is  0.00736  when  the  projectile  is  traveling  at  Mach  2  in  inviscid  helium  at  1  atm  and 
it  is  12%  higher  when  the  flow  is  viscous.  At  Mach  5,  the  drag  coefficient  increases  by  10% 
due  to  the  fluid  viscosity. 

Table  3:  Comparison  of  drag  force  per  unit  span  and  drag  coefficient  for  inviscid  and 
viscous  simulations  at  different  speeds.  Percents  are  with  respect  to  the  inviscid  case 
at  latm. 


Case 

Drag  Per  Unit  Span 

Drag  Coefficient 

Inviscid  1  atm 

Mach  2 

99N/m 

7.36x1 0’3 

Mach  5 

272  N/m 

3.24x1  O'3 

Mach  7 

681  N/m 

4. 13x1 0‘3 

Viscous  latm 

Mach  2 

1 1 1  N/m 

8.25x1 0'3 

(+12%) 

Mach  5 

302  N/m 

3.59x1 0’3 

(+10%) 

Mach  7 

770  N/m 

4.67x1 0’3 

(+13%) 

Inviscid  0.00  latm 

Mach  2 

0.0996  N/m 

7.40x1 06 

Mach  5 

0.277  N/m 

3.20x1  O'6 

Mach  7 

0.692  N/m 

4.20x1  O'6 

It  is  interesting  to  note  that  the  contribution  of  the  shear  drag  to  the  total  drag  on  the 
projectile  is  very  small  at  such  high  Reynolds  numbers  with  a  projectile  that  generates  only  a 
very  thin  wake.  If  the  channel  were  wider  and  the  shocks  were  not  to  bounce  back  on  the 
projectile,  the  viscous  effects  could  certainly  be  neglected.  However,  the  viscous  effects  are 
quite  significant  in  our  problem.  Because  of  the  presence  of  the  boundary  layer,  the  reflection 
on  the  projectile  walls  of  a  shock  is  radically  altered.  Hence  the  whole  shock  pattern  depends 
strongly  on  viscous  effects  at  the  walls,  and  so  does  the  net  aerodynamic  and  aerothermal 
forces  on  the  projectile. 
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Although  the  aerodynamic  heating  within  the  launch  tube  does  not  present  a  significant 
challenge  (due  to  the  low  pressure  and  density  of  the  gas)  the  heating  of  the  projectile  upon  exit 
of  the  launcher  into  the  dense  ground-level  atmosphere  is  an  issue  that  must  be  thoroughly 
examined  before  a  definitive  statement  on  the  feasibility  of  direct  launch  from  the  surface  of 
the  Earth  to  orbit  can  be  made.  For  example,  the  projectile  must  be  able  to  tolerate  extremely 
large  heat  fluxes  for  10-15  seconds  as  the  projectile  travels  at  around  7  km/s  through  the  dense 
lower  atmosphere.  Numerous  studies  have  indicated  that  flyout  velocities  of  around  6.5  km/s 
are  feasible  from  a  heat  transfer  perspective,  [1,2].  Furthermore,  these  references  indicate  that 
at  these  speeds,  “Sea  level  launch  poses  no  unusual  difficulty”,  [1].  Although  these  sources 
indicate  that  such  launch  speeds  are  feasible,  experimental  validation,  especially  using  ablative 
techniques,  has  not  been  performed.  Therefore,  immediate  work  has  begun  on  exploring 
appropriate  analytical,  numerical,  and  experimental  techniques  for  assessing  the  performance 
of  the  projectile  as  it  travels  through  the  dense  portion  of  the  atmosphere  at  extremely  high 
speeds  and  experiences  heat  fluxes  on  the  order  of  several  ten  kW/cm2,  and  a  significant  part  of 
the  Phase-2  effort  will  be  dedicated  to  optimizing  the  shape  of  the  projectile,  examining 
various  surface  ablator  materials  (such  as  plasma  sprayed  tungsten  which  has  recently  been 
perfected),  and  the  use  of  novel  multi-layer  aerodynamic  shells  to  withstand  the  high  heat 
loads.  Numerical  efforts  will  be  explored  to  compliment  the  work  currently  being  pursued  by 
the  University  of  Minnesota  and  the  University  of  New  Orleans. 


To  prepare  for  the  Phase  II  effort  in  understanding  what  types  of  protective  and/or  ablative 
coatings  are  necessary  to  protect  the  projectile,  it  is  first  essential  to  model  what  type  of  heat 
loads  the  projectile  will  actually  see  during  its  flight  through  the  atmosphere.  To  study  this 
issue,  a  model  was  developed  for  both  direct  orbital  insertion  and  flight  of  kinetic  energy 
projectiles  to  determine  what  heat  loads  are  expected  as  a  function  of  both  altitude  and  time. 
The  following  summarizes  the  model: 


■  As  from  the  nomenclature  used  in  Palmer&Dabiri2,  RB  is  the  base  radius  of  the 
conical  projectile  and  Rn  the  nosetip  radius. 

■  The  heat  flux  is  computed  from  Palmer&Dabiri  equation  (1),  and  is  assumed  to  be 
maximal  at  the  time  of  launch  since  the  ratio  of  atmospheric  densities  at  altitude  h 
and  launch  altitude  ph/po  are  then  maximal.  Thus,  for  a  launch  speed  V,  the 
maximum  heat  flux  is  taken  to  be 


(  V'8 

«3.4x1(T4(/?„)~0V8  . 

max  \  Po  ) 

■  The  launch  velocity,  V,  is  set  so  that  the  maximum  heat  flux  is  280  MW/m2. 

■  The  drag  coefficient,  Co,  is  assumed  to  be  only  a  function  of  the  cone  angle,  6C,  and 
the  ratio  of  nosetip  radius  to  base  radius,  according  to  Palmer&Dabiri  equation  (3), 
namely 


}_dQ 
A  dt 


=  2  sin2  + 


V’ 


K,J 


(j  -  2  sin2  9C 


2 

M.R.  Palmer  and  A.E.  Dabiri,  “Electromagnetic  space  launch:  a  re-evaluation  in  light  of  current 
technology  and  launch  needs  and  feasibility  of  a  near  term  demonstration”,  IEEE  Transactions  on 
Magnetics,  Vol.  25,  No.  1,  January  1989. 
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■  The  cone  angle,  0C,  is  4°. 

■  The  launch  angle  is  taken  to  be  the  one  for  which  the  projectile  has  no  vertical 
velocity  when  reaching  LEO  altitude  (180km  or  300km). 

■  Launch  takes  place  due  East  from  Kennedy  Space  Center. 

■  The  speed  required  for  circularization  at  LEO  1 80km  is  7,802m/s. 

■  The  speed  required  for  circularization  at  LEO  300km  is  7,732m/s. 

■  The  density  of  the  projectile  is  3000  kg/m3,  as  in  Palmer&Dabiri. 

■  The  mass  of  the  projectile  is  determined  from  its  volume  and  density. 


The  methodology  and  assumptions  for  these  computations  can  be  found  in  the  paper  by 
Uranga  et  al.3  In  particular,  the  circularization  speed  is  given  by 


V '  = 


gme 

RE  +  h 


go 


RE  +  h 


where  G  =  6.67  x  10“'‘m3/kg/s2  is  the  universal  gravitational  constant,  ME  =5.97x1 024kg  is 
the  mass  of  the  Earth,  RE  =  6.380  xl06m  is  the  radius  of  the  Earth,  h  is  the  orbital  altitude 
above  the  Earth’s  surface  in  meters,  and  g0  =  9.81  m/s2  is  the  acceleration  due  to  gravity  on  the 
Earth’s  surface. 


For  any  base  radius,  the  curves  in  Figure  25  and  Figure  26  have  a  maximum,  meaning 
that  there  is  a  ratio  of  radii,  Rb/Rn,  that  will  maximize  the  horizontal  velocity  at  LEO  altitude. 
In  other  words,  for  a  given  projectile  base  radius  RB  (so  for  a  given  projectile  volume  and 
hence  payload  size),  there  is  an  optimal  nosetip  radius  Rn  that  will  maximize  the  speed  at 
orbital  altitude  while  keeping  the  heat  flux  below  a  set  value  (280  MW/m2  here).  A  summary 
of  geometry  and  launch  parameters  for  launch  to  orbits  of  185  and  300  km  are  summarized  in 
Tables  2  and  3. 


3  A.  Uranga,  D.  Kirk,  H.  Gutierrez,  R.  Meinke,  and  K.  Barker,  “Rocket  performance  analysis  using 
electrodynamic  launch  assist”,  Proceedings  of  the  43rd  AIAA  Aerospace  Sciences  Meeting  and  Exhibit, 
Paper  No.  A1AA-2005-1449,  Reno,  Nevada,  Jan.  10-13,  2005. 
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Table  4:  Unpowered  projectile  to  LEO  (300  km). 


Rb 

(m) 

Mass 

(kg) 

Rn 

(m) 

Rb/Rn 

Launch  Mach 
Number  for 

280  MW/m2 

Launch 

Angle 

0 

Horizontal 
Velocity 
at  hLEO 
(m/s) 

0.01 

0.090 

0.0001 

100 

9.197 

65.9 

1,465 

0.000125 

80 

9.328 

64.5 

1,535 

0.0002 

50 

9.607 

62.2 

1,654 

0.00025 

40 

9.743 

61.3 

1,701 

0.0005 

20 

10.177 

60.7 

1,732 

0.001 

10 

10.631 

70.1 

1,260 

5.758 

100 

10.036 

47.4 

2,587 

80 

10.177 

46.6 

2,649 

50 

10.483 

45.1 

2,772 

40 

10.631 

44.4 

2,829 

20 

11.105 

42.8 

2,969 

10 

11.599 

42.9 

2,953 

8 

11.763 

43.9 

2,863 

5 

12.116 

49.6 

2,414 

4 

12.287 

55.7 

2,014 

89.963 

10.631 

42.3 

3,016 

80 

10.781 

41.6 

3,080 

50 

11.105 

40.2 

3,215 

40 

11.261 

39.6 

3,277 

20 

11.763 

37.9 

3,457 

10 

12.287 

36.9 

3,563 

5 

12.835 

38.6 

3,375 

2 

13.596 

59.4 

1,797 

40 


Figure  25:  Horizontal  projectile  velocity  when  reaching  LEO  (180  km)  versus  ratio  of  base  to 

nosetip  radii. 


Figure  26:  Horizontal  projectile  velocity  when  reaching  LEO  (300  km)  versus  ratio  of  base  to 

nosetip  radii. 
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The  model  also  calculates  optimized  projectile  geometry  for  a  given  base  radius  and 
desired  altitude  (or  range)  for  minimization  of  heat  flux.  This  is  in  contrast  to  most  models 
which  attempt  to  minimize  the  initial  muzzle  velocity.  In  other  words,  the  model  helps 
selecting  the  best  geometry  for  smallest  heat  load  to  the  projectile  for  a  given  mission. 

It  seems  feasible  to  utilize  experimental  methods  to  investigate  the  high  heat  fluxes 
originated  by  atmospheric  coasting  from  ground  by  making  use  of  electron  beam  welding 
techniques,  which  can  produce  surface  heat  fluxes  as  high  as  1 03- 1 06  kW/cm2,  [5].  Although 
the  actual  physical  situation  between  shear  layer  heating  during  atmospheric  flight  and  electron 
beam  penetration  to  achieve  the  heat  flux  is  quite  different,  exposing  various  nose  cone 
geometries  or  cooling  techniques  to  these  intense  heat  loads  will  nonetheless  provide  valuable 
guidance  for  the  design  of  the  actual  projectile. 

If  the  projectile  can  withstand  the  high  heat  loads  during  the  initial  ascent  (first  several 
seconds),  the  heat  flux  then  monotonically  decreases  to  much  more  tolerable  levels  (2-5 
kW/cm2).  Current  research  suggests  that,  “Lightweight  flexible  phenolics,  such  as  PhenCarb- 
28  at  28  Ib/fl3,  look  promising  for  [planetary  exploration]  missions  with  peak  heating  in  the 
range  from  about  1,000  to  2,000  Btu/ft2  -sec  (~1.1  to  2.3  kW/cm2).  A  Phencarb-32  or  -36  should 
have  more  optimal  performance  and  efficiency  for  still  higher  heating”,  [5].  The  existing  arc- 
heated  facility  -NASA  MSFC  PRL-  is  capable  of  producing  heat  fluxes  on  the  order  of  0.5 
kW/cm2,  so  heat  transfer  experiments  simulating  the  projectile  at  times  after  flight  through  the 
densest  portions  of  the  atmosphere  (after  the  first  1-5  seconds)  can  be  readily  accomplished. 

Another  option  includes  plasma  sprayed  high-temperature  refractory  metals  and 
ceramics.  For  example,  "Ultramet  has  developed  and  repeatedly  demonstrated  a  refractory 
ceramic  coating  material  capable  of  non-ablating,  long-term  operation  (minutes  to  tens  of 
minutes)  under  air  arcjet  conditions  at  heat  flux  and  enthalpy  levels,  simulating  leading  edge 
reentry  conditions,  of  350  Btu/ft2sec  (~0.4  kW/cm2)  and  12,000  Btu/lbm  respectively.  This 
unique  material  has  also  survived  heat  flux  levels  as  high  as  800  Btu/ft2 sec  (-1  kW/cm2), 
without  erosion,  under  lower  enthalpy  conditions.  This  patented  coating  material,  designated 
Ultra2000,  is  composed  of fine  alternating  layers  of  hafnium  carbide  (HfC)  and  silicon  carbide 
(SiC),  applied  by  chemical  vapor  deposition  (CVD),  to  a  total  thickness  of  just  0.005-0.010". 
Ultra2000  is  capable  of  long-term,  non-ablating  operation  at  surface  temperatures  as  high  as 
4000°F,  and  has  been  effectively  applied  to  and  tested  on  various  ceramic  matrix  composites, 
carbon-carbon  composites,  and  graphite  substrates”,  [5]. 

In  summary,  it  appears  that  if  the  projectile  can  survive  the  intense  heat  fluxes  for  the 
first  few  seconds  of  launch  (1-3  seconds  to  arrive  at  16,000  km,  where  the  atmospheric  density 
has  already  dropped  significantly),  proven  thermal  protection  systems  are  then  available  to 
protect  the  projectile  for  the  remainder  of  the  flight.  The  Phase-2  contract  will  explore  each  of 
these  options,  and  use  appropriate  teaming  (Ultramet,  Plasma  Processes,  Inc.,  etc.)  to  develop 
and  test  an  appropriate  thermal  protection  system. 
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7.  Combined  Simulation  of  Electromagnetic  and  Aerodynamic  Forces 

An  integrated  simulation  tool  that  combines  the  electromagnetic  and  aerodynamic  forces 
acting  on  the  projectile  is  currently  being  developed  to  calculate  the  multi-DOF  projectile 
trajectory  during  acceleration.  The  simulation  consists  of  the  following  sub-modules: 

1.  Computation  of  the  instantaneous  electromagnetic  forces  and  torques  induced  in  the 
projectile  by  the  currently  active  track  coils  (on  Figure  1,  Tri  and  7>/).  In  the  example’s  rear 
section,  these  depend  on  the  currents  /*/,  />/,  the  local  axial  position  x,  and  the  geometric 
relationships  between  track  and  projectile. 

This  module  consists  of  the  following  sub-modules: 

1.1  Computation  of  the  instantaneous  mutual  inductance  between  the  projectile  coils  and 
the  corresponding  track  coils.  In  the  example  this  means  Tri  with  Pr  and  7>/  with  Pp.  Inputs  to 
this  sub-module  are: 

•  The  coils’  radii  (Rp  and  R,  for  projectile  and  track,  respectively). 

•  The  <x,y,z>  vector  of  relative  position  between  coil  centers  (x=  distance  between 
coils’  centers  along  the  launch  axis,  y=  horizontal  distance,  z=  vertical  distance). 

•  The  pitch  ( Q )  and  yaw  (9Z)  angles  of  both  projectile’s  coils  relative  to  the  local 
coordinate  system  X fixed  on  the  corresponding  track  coil. 

Single  turn  coils  represented  by  infinitely  thin  filaments  will  be  considered  first.  Later, 
approximate  3D  coil  geometry  (coil  thickness  and  depth)  will  be  included. 

1.2  Computation  of  the  instantaneous  electromagnetic  propulsion,  lift  and  guidance 

forces  acting  on  the  projectile.  This  sub-module  consists  of: 

•  Evaluation  of  the  gradient  of  the  mutual  inductance  calculated  in  Module  1.1  for  the 
current  <x,y,z>  vector  of  relative  position  between  coil  centers. 

•  Computation  of  the  <Fx,Fy,Fz>  vector  of  electromagnetic  forces  acting  on  both  front 
and  rear  projectile  coils.  This  is  based  on  the  mutual  inductance  gradient  calculated 
above  and  the  corresponding  instantaneous  currents  (/*/  with  ipR,  ipj  with  ipp). 

1.3  Computation  of  electromagnetic  torques  acting  on  the  projectile.  The  forces 
calculated  on  Module  1.1  will  be  used  to  calculate  the  electromagnetic  torques  created  by  the 
propulsion,  guidance  and  lift  forces. 

2.  Computation  of  the  instantaneous  aerodynamic  forces  and  torques  induced  in  the  projectile 
during  launch.  This  module  consists  of  the  following  sub-modules: 

2.1  Computation  of  the  3D  aerodynamic  forces  created  by  the  pressure  distributions 
acting  on  the  projectile’s  surface,  as  a  result  of  the  reflection  of  shock  waves  within  the  launch 
tube.  A  simplified  model  of  the  pressure  distributions,  validated  by  computational  fluid 
dynamics  state-of-the-art  software,  will  be  used. 

2.2  Computation  of  aerodynamic  torques  acting  on  the  projectile  as  a  result  of  the  forces 
calculated  in  Module  2.1. 
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3.  Computation  of  the  6-DOF  trajectory  of  the  projectile  during  launch.  This  module  consists  of: 

•  The  electromagnetic  and  aerodynamic  forces  and  torques  calculated  by  Modules  1.2, 
1.3,  2.1  and  2.2  are  used  to  estimate  the  multi-DOF  projectile  trajectory  (based  on 
integration  of  Newton’s  equations  of  motion)  in  a  sub-section  of  the  track  such  as  the 
shaded  area  on  Figure  9. 

•  An  algorithm  that  connects  the  motion  in  the  present  sub-section  of  the  track  with  the 
next  sub-section. 

The  integrated  simulation  will  be  a  novel  and  valuable  tool.  It  will  be  used: 

•  To  prove  the  feasibility  and  design  requirements  of  the  CFSM  launcher  concept. 

•  To  develop  design  optimization  of  the  launcher  geometric  parameters. 

•  To  provide  insight  in  the  interaction  of  electromagnetic  and  aerodynamic  effects  during 
launch. 

During  Phase-2,  the  launcher  simulation  will  be  integrated  with  mathematical  models  of  the 
energy  sources,  switching  devices  and  other  power  electronics.  This  will  allow: 

•  To  estimate  the  performance  of  the  CFSM  concept  in  a  full-scale  design. 

•  To  optimize  parameters  of  the  power  conversion  system. 

•  To  study  the  need  and/or  feasibility  of  a  real-time  feedback  system  to  optimize  the 
performance  of  the  CFSM  system. 


44 


8.  Experimental  Qualification  of  the  CFSM  Concept 

The  constant-flux  principle  is  based  on  the  static  relationships  of  current-field-flux  given 
by  the  geometric  relationships  between  track  coils  and  projectile  coils.  Therefore  the  principle 
can  be  tested  under  semi-static  conditions  (low  speed).  Figure  5  shows  a  section  of  a  test  set-up 
used  to  demonstrate  the  constant-flux  principle. 


Figure  27:  Linear  test  rig  for  qualification  tests  of  the  CFSM  system. 

The  linear  test  rig  represents  a  section  of  the  launch  tube  with  a  few  track  coils 
implemented.  A  projectile  mock-up  will  be  moved  along  the  tube  to  trigger  the  sequence  of 
current  pulses  f2(t),  fi(t).  The  mock-up  projectile  is  rigged  with  a  flux  sensing  coil  -  any 
flux  change  experienced  by  the  test  coil  would  induce  a  voltage  that  can  be  externally 
measured.  Integrating  the  signals  from  the  flux  sensing  loops  measures  the  total  flux  change 
seen  by  the  test  coil.  The  projectile  mock-up  can  be  moved  to  any  position  relative  to  the 
pulsed  track  coils,  and  the  current  waveforms  calculated  by  solving  (11)  subject  to  (12)  can  be 
generated  by  computer-controlled  current  sources.  Applying  the  calculated  pulse  shapes  to  the 
mock-up  track  coils  and  measuring  the  induced  voltage  simulates  the  flux  seen  by  the  coil  on 
the  moving  projectile.  Flux  measurements  as  a  function  of  time  will  be  performed  for  all 
positions  between  two  track  coils.  For  a  correctly  shaped  current  pulse,  the  measured  voltage 
and  therefore  the  flux  should  be  constant  at  all  test  positions. 

The  tests  will  show  what  level  of  accuracy  in  pulse-shaping  is  needed  to  experience 
acceptable  flux  changes  in  the  actual  system,  or  otherwise,  which  flux  changes  will  be  seen  by 
the  moving  projectile  coils  for  a  given  inaccuracy  in  the  actual  pulse  shape.  These  tests  and  the 
corresponding  calculations  will  help  optimize  the  CFSM  system,  and  estimate  the  induced 
voltages  (and  therefore  parasitic  currents)  induced  in  the  moving  projectile  as  a  result  of  non¬ 
ideal  pulse  shaping.  While  it  is  in  principle  possible  to  adjust  the  pulse  shapes  to  any  level  of 
accuracy,  system  complexity,  robustness  and  cost  suggests  the  use  of  the  simplest  pulse  power 
system  that  will  not  cause  quenching  and  overheating  of  the  projectile  coils.  Various  coil 
configurations  for  track  and  projectile  can  be  tested  and  compared  with  calculations.  The 
linear  test  rig  could  also  be  used  to  measure  the  propulsion  and  levitation  forces,  using  active 
coils  in  the  projectile  mock  up.  Since  the  mock-up  would  only  move  over  short  distances, 
external  connections  for  the  coils  would  be  easy  to  implement. 
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8.1  Superconductor  Qualification  Tests 
Quench  Stability 

A  prototype  coil  consisting  of  a  few  turns  of  superconductor  would  be  excited  to  the 
required  operational  current.  A  second  coil  in  close  proximity  would  be  used  to  simulate  the 
flux  changes  that  are  expected  in  actual  operation  of  these  coils.  Inducing  flux  changes  of 
various  levels  into  the  prototype  projectile  coil,  its  quench  stability  can  be  tested.  The  flux 
changes  used  in  these  tests  can  be  stronger  than  expected  in  the  actual  system  to  guarantee  a 
sufficient  operational  margin  during  actual  operation. 

Mechanical  Stability 

During  launch,  the  superconductor  in  the  projectile  coil  will  experience  significant 
mechanical  forces  due  to  inertial  acceleration.  Superconductors  (either  LTS  or  HTS)  have 
never  been  tested  under  such  conditions.  The  inertial  acceleration  is  equivalent  to  compression 
of  the  conductor,  i.e.,  the  mechanical  launch  conditions  can  be  simulated  by  compressing  the 
conductor.  To  test  stability  to  mechanical  stress,  the  conductor  would  be  placed  in  a  cryostat 
with  a  background  field  of  appropriate  strength  that  simulates  the  self-field  of  the  projectile 
coil  and  the  field  due  to  the  track  pulses.  The  conductor  is  then  loaded  to  the  appropriate 
mechanical  pressure  and  its  critical  current  measured.  The  conductor  and  its  enclosure  should 
pass  these  tests  without  quench  with  at  a  significant  safety  margin.  These  tests  can  be 
performed  with  different  superconductors  to  identify  the  best  choice. 

8.2  Study  of  Nose  Cone  Ablation 

As  already  presented  in  Section  6,  the  coasting  flight  of  the  projectile  through  the 
atmosphere  remains  a  major  issue  and  requires  special  testing.  High  intensity  particle  beams,  as 
used  on  electron  beam  welding,  produce  sufficient  heat  flow  densities  to  simulate  the  effects  of 
low-altitude  atmospheric  coasting  on  the  projectile. 

It  seems  feasible  to  use  experimental  methods  to  investigate  the  high  heat  fluxes  originated 
by  atmospheric  coasting  from  ground  by  making  use  of  electron  beam  welding  techniques, 
which  can  produce  surface  heat  fluxes  as  high  as  103-106  kW/cm2  [3].  Electron  beams  used  for 
welding  are  of  rather  lower  energy,  and  the  particles  are  absorbed  in  thin  layer  of  the  irradiated 
object.  However,  the  drag  form  the  atmosphere,  which  might  quasi  peel  off  the  heated  layer 
would  be  difficult  or  impossible  to  simulate  in  such  experiments.  Spinning  the  irradiated  object 
might  be  a  possibility  to  simulate  such  a  peel-off  effect  from  the  atmosphere.  Although  the 
actual  physical  situation  between  shear  layer  heating  during  atmospheric  flight  and  electron 
beam  penetration  to  achieve  the  heat  flux  is  quite  different,  exposing  various  nose  cone 
geometries  or  cooling  techniques  to  these  intense  heat  loads  will  nonetheless  provide  valuable 
guidance  for  the  design  of  the  actual  projectile. 
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9.  System  Integration 

A  conceptual  design  of  a  complete  electromagnetic  launch  system  for  affordable,  agile 
access  to  space  consisting  of  a  approximately  500-m  long  track  with  a  inclination  of  about  45 
degrees,  the  required  propulsion  system,  the  reduced-pressure  helium  atmosphere  in  the  launch 
tube  has  been  performed  during  Phase- 1  and  is  described  in  previous  and  the  current  Status 
reports. 

The  tests  described  in  section  4  of  this  report  will  enable  us  to  perform  a  more  detailed 
design  of  the  system  during  Phase-2.  This  design  will  also  include  the  required  pulse  power 
system  and  the  energy  storage  system,  and  a  concept  for  the  low  pressure  launch  tube.  Since 
large  forces  will  be  experienced  by  the  track  coils  of  the  CFSM  system,  a  concept  for 
supporting  these  will  be  presented.  The  track  support  concept  has  to  include  requirements  for 
the  alignment  of  the  launch  tube.  Special  attention  has  to  be  given  to  the  fact  that  the  launch 
tube  has  to  consist  of  a  non-conductive  material,  most  likely  a  fiber-reinforced  composite. 

9.1  Scalability 

The  proposed  launch  system  concept  has  great  potential  for  scalability.  Although  our  design 
studies  focused  on  a  projectile  mass  of  about  100  kg,  consisting  of  a  launch  vehicle  with  a 
mass  of  5-10  kg  and  a  roughly  10  times  larger  support  structure  with  superconducting  coils  in 
their  cryostat  and  the  nose  and  rear  cones  (material  of  cones  is  TBD,  see  section  4.1  and  4.2)  of 
the  projectile,  it  seems  possible  to  design  similar  systems  with  different  parameters.  Based  on 
the  work  performed  so  far,  it  seems  possible  to  increase  the  projectile  mass  within  realistic 
limits,  in  particular,  if  the  requirements  for  acceleration  (currently  5000  gee)  and/or  the 
required  launch  velocity  (currently  7000  m/sec)  are  reduced.  Such  variations  of  design 
parameters  should  open  up  possibilities  for  various  other  applications  of  the  system. 
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and  the  advantage  of  helium  versus  air  is  reduced. 
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Optimization  of  Acceleration  Profile 
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1  200  400  600 

Position  on  Track  [m] 


The  differential  equation  describing  the  acceleration  profile  can  be  solved  analytically  for  velocity  and  position  of  the  sled  as  a  function  of  time  t. 
^start 's  acceleration  at  t  =  0  and  m0  is  the  slope  of  the  acceleration  profile  (see  avove).  The  following  solution  is  obtained. 
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I  300  400  500 

Track  Position  [m] 


Determine  sled  velocity  as  a  function  of  track  position 


Track  Position  [m] 


The  required  power  as  a 
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Define  acceleration  profile  bstart  :=  5550-g  bend  :=  4100-g  Corresponding  length  for  acceleration  change:  As  :=  500-m 


a 
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Estimate  total  force  including  , 

magnetic  drag  (  ~5  %)  bAcc(s)  •—  1 .05- F(0t (s)  FAccC^’in)  =  6  x  10  kN  Fj\cc(S>a)  —  4.37  x  10  kN 


00 

00 


Total  Propulsion  forces  in  Helium 


The  diagram  shows  the  required  force  (blue  solid 
curve)  for  helium  with  the  same  acceleration  profile  as 
for  air.  Since  the  drag  in  the  helium  atmosphere  is 
smaller  than  in  air,  the  required  force  falls  off  towards 
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The  required  power  as  a  function  of  track  position  for  the  new  (optimized)  acceleration  profiles  is  shown  in  the  following  diagram: 


00 
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General  definitions: 


Drag  coefficient  versus  Mach  number 


coefficient  versus  Mach  number  can  be  used  for  other  values  of  the  drag  coefficient.  Under  these  assumptions  the  velocity  dependence  of  any 
drag  coefficient  can  be  estimated  by  multiplying  the  coefficient  with  the  function  draglncr(v,  T,  y).  The  argument  v  is  the  velocity  in  m/sec,  T  is 
absolute  temperature  in  K  and  y  is  the  heat  capacity  ratio  Cp/Cv  ~  1 .4  (air)  and  1 .67  (helium).  R  =  287  J/kg/K  is  the  ideal  gas  constant. 


